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ANALYSIS OF HORIZONTAL- TAIL LOADS IN PITCHING MANEUVERS 
ON A FLEXIBLE SWEPT-WING JET BOMBER 
By William S. Aiken, Jr. 

SUMMARY 


Horizontal-tail loads measured by means of strain gages in pitching 
maneuvers are analyzed to determine wing-fuselage aerodynamic-center posi- 
tion, zero-lift pitching -moment coefficient, airplane pitching moment of 
inertia, and radius of gyration. A similar analysis is made of the time- 
history data for the elevator angles and the results were found to agree 
with those from the tail-load analysis. The flight-determined values of 
aerodynamic-center position for rigid conditions and the zero-lift 
pitching-moment coefficients were in some disagreement with the wind- 
tunnel data over the Mach number range of the tests (0.42 to 0.8l). The 
pitching moment of inertia determined from the flight data for rigid-wing 
conditions agreed with calculations based on ground tests. The effective 
pitching moment of inertia computed from theoretical consideration for 
flexible flight conditions was in disagreement with flight data. Details 
of the analysis procedures and least-squares methods used are given. 


INTRODUCTION 


The calculation of airplane design tail loads and stability charac- 
teristics requires reliable estimates of the wing-fuselage pitching-moment 
characteristics. The use of highly swept flexible wings combined with 
other flexible airplane components introduces additional factors which 
must be considered in tail-load design analysis procedures. Investiga- 
tions by the National Advisory Committee for Aeronautics of a large 
flexible swept-wing jet bomber which included measurements of horizontal- 
tail loads permitted the analysis of data from which comparisons could 
be made between wind-tunnel measurements of wing -fuselage aerodynamic - 
center positions and zero-lift pitching-moment coefficients and values 
of these parameters as derived from flight data. 

The analysis of flight data in the present report is, to a large 
extent, based on analyses and information contained in references 1 
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and 2 for wing deflections, reference 3 for horizontal-tail parameters, 
reference 4 for airplane lift-curve slopes and angles of zero lift, and 
reference 5 for wing centers of pressure. The methods used to analyze 
the flight data and to convert measured pitching-moment parameters to 
equivalent rigid conditions for comparison with wind-tunnel data are 
described in detail. Comparisons are given between flight and wind- 
tunnel results for aerodynamic -center position and zero-lift pitching- 
moment coefficients and between flight and calculated values of moments 
of inertia and radii of gyration for both flexible and rigid conditions. 
Although no direct comparisons are made with present tail-load design 
computation methods, the theoretical methods which were used in the 
flight-data analysis contain the essential elements of design procedures 
for flexible aircraft and, thus, provide an indirect check on their 
adequacy. 


SyMBOLS 


A,B,C,D coefficients of equations (A7) and (Bl) used to obtain air- 
plane pitching-moment parameters 


A corr the A coefficient of equation (Bl) corrected for zero 

shift 



horizontal-tail lift-curve slope per degree 



horizontal-tail lift-curve slope per degree for flexible 
fuselage conditions defined by the expression 





horizontal-tail lift -curve slope per degree with root eleva- 
tor angle 


C, 


tail pitching -moment coefficient due to elevator deflection 


C, 


•mo 


zero-lift wing- fuselage pitching-moment coefficient 


C, 




zero-lift wing-fuselage pitching-moment coefficient com- 
puted directly from measured zero-lift tail load 


C Mqj 


corrected zero-lift wing-fuselage pitching-moment coeffi- 
cient using analysis method I 
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Cmojj zero-lift wing-fuselage pitching-moment coefficient using 

analysis method II 

Cjy. airplane normal-force coefficient corrected for pitching- 

acceleration tail load 

Iy airplane pitching moment of inertia, slug -ft ^ 

Iyj. effective airplane pitching moment of inertia defined by 

equation (AlO) , slug-ft2 

aerodynamic tail load, lb 

Lt aerodynamic tail load plus component of tail aerodynamic 

pitching moment defined by equation (a 6), lb 

Iy aerodynamic wing load, lb 

Ly-.- aerodynamic wing load per unit pitching acceleration, 

9 

lb 

radian/sec^ 

La aerodynamic wing load due to pitching velocity (eq. (A3)), 

lb 

Lq aerodynamic wing load due to pitching acceleration (eq. (A2)), 

lb 

M Mach number 

M ac pitching moment about wing-fuselage aerodynamic center 

N number of equations in least-squares solutions 

S wing area, sq ft 

S-t horizontal-tail area, sq ft 

Tt tail aerodynamic torque, in-lb 

V true airspeed 

W airplane weight, lb 

horizontal-tail weight, lb 
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location of airplane center of gravity, percent c 

location of rigid-wing — fuselage additional-load aerodynamic 
center, percent c 

location of flexible-wing — fuselage additional- load aerody- 
namic center, percent c 

location of center of pressure of aerodynamic wing load due 
to pitching acceleration, percent c 

zero shift in measured aerodynamic tail loads (from ref. 3)> 
lb 

zero shift in measured aerodynamic tail torque, in-lb 

faired tail-on airplane lift -curve slope per degree (from 
ref. k) 

wing mean aerodynamic chord, in. 
tail mean aerodynamic chord, in. 

distance from wing-fuselage aerodynamic center to airplane 
center of gravity, positive with center of gravity forward 
of the aerodynamic center, in. 

center of pressure of wing load due to pitching acceleration 
in. 

center of pressure of wing load due to pitching velocity, in 
acceleration of gravity, ft/sec^ 
tail incidence, deg 

wing incidence, deg 

airplane radius of gyration in pitch, ft 
effective airplane radius of gyration in pitch, ft 

horizontal-tail length, distance from airplane center of 
gravity to quarter-chord of horizontal-tail mean aerody- 
namic chord, in. 
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rigid-wing— fuselage lift-curve slope per degree 
normal load factor at airplane center of gravity- 
normal load factor at horizontal tail 
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m R 

n 

“t 

q. 

w 

X, 


AX, 


ac 




°t 

“v 


Xdj 


e i 


e 2 

0 

0 


I 


dynamic pressure, 


lb/sq ft 


weighting factors 

distance from wing-fuselage aerodynamic center to quarter- 
chord of horizontal-tail mean aerodynamic chord, in. 

difference between rigid- and flexible -wing— fuselage 
aerodynamic— center positions, defined by equation (6), 
percent c 

2 

difference between theoretical and measured ky values, 
ft 2 


tail angle of attack, deg 


wing angle of attack, deg 

wing angle of zero lift (from ref. 4), deg 

root elevator angle, positive down, deg 
root elevator angle at zero load factor, deg 

errors in fit or measurements; subscripts to e denote 
error associated with quantity indicated 

2 

error in calculated ky^ 

p 

error in measured ky^ 
pitching velocity, radian/sec 

measured pitching acceleration corrected for instrument 
response characteristics, radian/sec^ 

indicates summation 
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de/da effective downwash factor 

Matrix notation: 

column matrix 

rectangular matrix 

transpose of rectangular matrix 

inverse of square matrix 

Bars over symbols indicate average values. 

APPARATUS AND TESTS 



Airplane 

The airplane used for this investigation was a six-engine, jet- 
propelled medium bomber. A photograph of the test airplane is shown as 
figure 1, and the airplane and horizontal-tail characteristics and dimen- 
sions are given in table I. 


Instrument at ion 

The data used for analysis in the present paper were obtained from 
standard NACA recording instruments and from strain gages mounted on the 
right and left sides of the horizontal tail. 

Normal accelerations were measured by two air-damped accelerometers, 
one near the center of gravity and one at the horizontal tail. Angular 
accelerations in pitch were measured by a rate-gyro type, electrically 
differentiating, magnetically damped turnmeter. Airspeed and altitude 
measurements were made with an NACA pitot— static head mounted on a boom 
approximately 1 maximum fuselage diameter ahead of the original nose. 

Electrical wire-resistance strain gages (Type A-6) with low tempera- 
ture correction factors were used to measure the root shears, bending 
moments, and torques at stations on the right and left sides of the tail. 
The gages were installed as four-active-arm bridges on the web and flanges 
of the main spars (50 percent chord) and on the upper and lower skin sur- 
faces near the leading edges of the horizontal tail. 
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The strain-gage bridge installation was calibrated according to the 
method detailed in reference 6. The bridges were then combined electrically 
so that, except for secondary carryover effects, a combined shear, moment, 
or torque bridge responded primarily to the shear, moment, or torque for 
the side of the tail on which the load was being measured. Final cali- 
bration equations using combined bridge outputs included carryover-effect 
corrections . 

The combined strain-gage outputs were recorded on an 18-channel 
oscillograph with individual galvanometer responses flat to 60 cps. All 
data were evaluated by using nondimens ional deflections as 


_ Flight deflection - Ground zero deflection 
P Calibrate signal deflection 

The sensitivity of each combined bridge was generally recorded prior to 
entering a maneuver through the use of a calibrate signal. With this 
system of data reduction, changes in battery voltage had no effect on the 
measurement of loads. In addition, galvanometer zeros with strain— gage 
power off were recorded to compensate for any mechanical shifts in the 
galvanometer zero position due to temperature effects in the recorder 
and any thermal electromotive— force effects in the strain— gage circuits. 

Aerodynamic tail loads on the horizontal tail were obtained from 
the structural loads (measured by the strain-gage bridges) and the known 
tail weight and normal load factor from the equation 

^aero " "^struct ” °t^t 

The aerodynamic bending moments and torques were obtained in a similar 
manner . 

The recorded data for all instruments were synchronized at 0.1-second 
intervals by means of a common timing circuit. All instruments were damped 
to about O.67 of critical damping. A summary of pertinent quantities meas- 
ured, instrument locations, and accuracies are given in the following 
table: 
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Quantity 

Location 

Instrument 

range 

Instru- 

ment 

accuracy 

Normal acceleration, 
g units ...... 

3^.2 percent c 

0 to 2 

0.005 

Normal acceleration, 
g units 

47.8 percent horizontal- 
tail root chord 

-2 to 6 

0.02 

Pitching acceleration, 

radians/sec 2 . . . 

25 percent c 

± 0.50 

0.01 

Dynamic pressure, 
lb/sq ft 

140 in. ahead of 
original nose 

0 to 800 

1 

Static pressure, 
lb/sq ft 

132 in. ahead of 
original nose 

0 to 2,200 

2 

Tail shear, per 

side, lb 

Root of tail 

± 25,000 

60 

Tail torque, per 
side, in- lb .... 

Root of tail 

± 2 , 000,000 

4,000 


Tests 

All tests were made with the airplane in the clean condition. The 

flight data evaluated in this report were taken from 68 push-down pull- 

up maneuvers (the same maneuvers used in refs. 3 and 4 ) made at alti- 
tudes of 20,000, 25,000, 30,000, and 35,000 feet and an overall Mach num- 
ber range from 0.427 to 0.812. The tests were made at normal and forward 
cent er-of -gravity positions and airplane weights ranging from 104,000 
to 127,000 pounds. Table II is a summary of the flight conditions for 
these runs. In the table are listed the flight and run numbers, average 
Mach number, average dynamic pressure, test altitude, weight, and center- 
of -gravity position. The range of Mach number and dynamic -pres sure changes 
during any test are also indicated. It might be noted that the center-of- 
gravity listings in table II differ slightly from those given for the same 
maneuvers in references 3 and 4. The airplane centers of gravity have been 
corrected for the effect of airplane attitude on the fuel-tank centers of 
gravity for the three large unbaffled fuselage tanks. 
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METHODS AND RESULTS 


The airplane pitching -moment parameters, that is, wing-fuselage 
aerodynamic-center position and zero-lift pitching -moment coefficient, 
and the airplane effective moment of inertia may be evaluated from flight 
measurements of the airplane motions and the horizontal-tail load. In 
appendix A, use is made of the airplane pitching-moment equation to set 
up two methods of analysis amenable to least-squares treatment. The 
methods are: 

Method I. A procedure in which direct tail-load measurements are 
used. 

Method II. A procedure in which elevator -angle measurements are 
used. 

Pitching-moment parameters obtained by the use of either method include 
quasi-static wing-flexibility effects but the equations do not allow for 
dynamic wing, tail, or fuselhge frequency-response effects since the bulk 
of the data presented was obtained without excitation of the major air- 
plane components. 

The following sections present the determination of the pitching- 
moment parameters from flight time-history data and comparisons with 
available wind-tunnel and mass-distribution data. The method used for 
extrapolating or correcting the measured aerodynamic -center positions 
to rigid-wing conditions is given in detail. The theoretical relation- 
ship existing between the measured or effective moment of inertia and the 
actual or rigid-airplane moment of inertia is also described. 


Basic Data 

The least-squares data-reduction procedures as used for methods I 
and II are given in appendix B. These procedures are used for the evalua- 
tion of the required tail-load and elevator -angle coefficients for each of 
the 68 maneuvers studied. These tail-load and elevator -angle coefficients 
in turn are used for the evaluation of the airplane pitching-moment param- 
eters. Both methods are illustrated by use of the time-history data for 
n, 0, Ltj, and & shown in figures 2, 3, and k for an example 

maneuver (flight 12, run 27). The calculated tail-load and elevator- 
angle time histories from equations (b 4) and (b 6) are illustrated in fig- 
ures 3 and if, respectively. 

It is demonstrated in appendix B that for method I a simplified form 
of equation (B2) which omits the 01^/v term could be used for the 
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determination of the zero-lift tail load, the tail load per g, and the 
tail load per unit pitching acceleration. The A, B, and C coeffi- 
cients of equation (B2) , along with their standard errors and errors of 
fit, for each of the 68 maneuvers analyzed are given in table III. 

The coefficients from the least -squares analysis of the elevator- 
angle data (method II), along with their standard errors and errors of 
fit, for each of the 68 maneuvers analyzed are given in table IV. 


The B coefficients of method I and the 88 /8n coefficients of 
method II, given in tables III and IV, respectively, as deduced from 
flight time-history data now permit the determination of the wing-fuselage 
aerodynamic -center location. This section illustrates the methods used to 
extract the aerodynamic-center data and to extrapolate the data as meas- 
ured for flexible-wing conditions to rigid-wing conditions. In appendix A, 
equations (A8) and (A25) show the relationship between the aerodynamic- 
center position d (the distance between airplane center of gravity and 
wing-fuselage aerodynamic-center location) and the measured coefficients. 

A general equation expressing the aerodynamic-center position in terms of 
its location on the wing mean aerodynamic chord is 


and is used to correlate the data obtained at various center- of -gravity 
locations . 

Aerodynamic -center position using method I .- Equation (A8) of appen- 
dix A gives the aerodynamic-center location d as 


Inserting numerical values for the example maneuver of appendix B with 
B = 392 (from eq. (b4)), l t = -552 inches, and W = 110,300 pounds 
in equation (2) gives 


The aerodynamic center in terms of the mean aerodynamic chord, using 
equation (l) and the cent er-of -gravity position of 22.9 percent, thus 
becomes 


Aerodynamic-Center Position 



( 1 ) 


d = -1.97 inches 
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fee) flex = 22 *9 + = 21 ‘ 6 P ercent 

The error in (X ac ) flex may he determined by use of the standard error 
in the B coefficient as 


% 


_ e B Z t 100 
c ~ W - B c 


(3) 


Using the standard error of B from table III gives 

358(-552) 


%c = ± 


(109,910) (1.559) 


= ±1.2 percent 


The aerodynamic-center positions and the associated standard errors 
for each of the 68 maneuvers are given in table V. 

Aerodynamic-center position using method II ,- Equation (A25) of 
appendix A gives the aerodynamic -center location d as 



”d6 da t + di t h de + L de\ W 

dn db dn-t y2 da \ da/a^qS 

+ ^C c 

an S t 

qS t “ (Sjf 

S6 da t + di t z t de + L de\ W 

an d& dn-t v 2 da \ da/apqS 



( \ ua+- de 

and Cmg were obtained from reference 3 and the airplane lift-curve 
slope ap was obtained from reference 4. The parameter (tail- 

incidence change due to fuselage bending under inertia loads) was obtained 
from equation (7) of reference 3. The remaining parameters required are 
g, the acceleration of gravity, V, the true airspeed, and q, the dynamic 
pressure . The quantity (36 /dn is the coefficient of equation (A19) 
associated with n and is given in table IV. 

For the example maneuver, substitution of numerical values into 
equation (4) results in 


d = -3-22 inches 
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The aerodynamic center by method II for the example maneuver is from 
equation (l): 


( X ac) flex 22 * 9 + (^ 559 , 


= 20.8 percent 


The error in (X a c)flex M y be computed from the standard error of the 
coefficient given in table IV as 



H( 

; c s) f I? + 

€n k 100 
e d&/dn 5 

<l s t ~ ( Clt 0 f 

& + lit + H de + L _ de\ W 

dn d 8 dn-t y-2 da \ da,/apqS 


(5) 


Using the standard error of 86 / 8 n from table IV gives 


[(-552) (0.0559) (0.420) + (-0.0086)105] (±0.258) 

ex — ■ — ■ — — = ±0.9 percent 

ac 1.559(2.562) 

The aerodynamic-center positions and their associated standard errors 
for all 68 maneuvers are given in table V where they may be compared with 
the values determined by using method I. 

Extrapolation to rigid-wing conditions .- The aerodynamic -center- 
position data in table V are for flexible-wing conditions. If the effects 
of wing flexibility are known, the flight measurements may be extrapolated 
to rigid-wing conditions and the variation of aerodynamic-center position 
with Mach number established. Data are available in reference 5 which 
may be used for this extrapolation. The forward shift of the wing 
aerodynamic-center position as a function of the flexibility parameter qmp 

is shown in figure 5> as determined from the theoretical curve (for an 
average value of W = 110,000 pounds) of figure 4(c) in reference 5 by con- 
version of the root center-of -pressure variation with qm-p to percent mean 

aerodynamic chord. Since the wing additional-load center-of -pressure 
data in figure 4(c) of reference 5 ^ere determined from wing-root aero- 
dynamic torques and the airplane center-line shear nW - i , the 

rigid-wing — fuselage aerodynamic -center position may be determined from 
the equation 
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(^ac)r " (^ac)flex + ^ac 


(6) 


When equation (6) is used it is assumed that the wing-fuselage aerodynamic- 
center position differs from the wing aerodynamic-center position only 
by a constant for any given Mach number and that changes in dynamic pres- 
sure do not p,ffect the fuselage contribution to the total wing-fuselage 
aerodynamic -center position. 

The values of qm R for use in the determination of from fig- 

ure 5 are listed in table V and were obtained from the dynamic pressures 
given in table II and the rigid-wing— fuselage lift-curve slopes m R 

given in reference 4. The (Xa C ) r values obtained from equation (6) for 
the method I and method II data and the AX^ values from figure 5 are 
given in table V. The values of (^ac)r methods are plotted 

in figure 6 and are identified by method. 

It will be noted that the errors associated with the aerodynamic 

centers ev of table V are not constant. Use was made of these errors 
A ac 

to define weighting factors to obtain weighted average values of aero- 
dynamic center at the group Mach numbers indicated in table V. The 
weighting factor is defined as 


w = 


; *ac 


( 7 ) 


and the weighted average aerodynamic-center position is defined by the 
equation 

M, ■ 1 


( 8 ) 


z 


w 


The standard error of the weighted average is given by the equation 


£ *ac 


^ w ( x ac)r 2 - (*ac)r w ( x ac)i 


" z 


(9) 


w 


The last column of table V gives the weighted average values of the 
aerodynamic-center position and the standard errors computed using equa- 
tions (8) and (9) for the Mach number groups used. These weighted average 


» 
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aerodynamic-center values are also plotted as circles in figure 7 with 
a line faired through the data to indicate a reasonable variation with 
Mach number, the standard error of each point (shown as the vertical 
lines within the Symbols) being considered. 

Comparison with wind-tunnel data .- Wind-tunnel model data corrected 
for model flexibility effects are given in reference 7 on pages L-124, 
L-126, and L-128. From these data the wing-fuselage aerodynamic-center 
data shown in figure 7 (the diamond- shaped symbols) were obtained for the 
available Mach numbers and may be compared with the flight measurements. 

Aerodynamic-center positions at various altitudes .- The wing -fuselage 
aerodynamic -center position as affected by wing flexibility may be calcu- 
lated for various altitudes by the use of equation (6), the AX ac data 

of figure and the faired curve of figure 7- The results of these cal- 
culations are shown in figure 8 and are considered to be the best esti- 
mates of aerodynamic -center position that can be made from the flight 
data. The results are limited to a low Mach number of 0.40, a high qm^ 

of and an airplane weight range from 110, 000 to 130,000 pounds since 

these conditions represent the limits within which the flight measure- 
ments were obtained. 


Zero-Lift Pitching -Moment Coefficient 

In the following section the determination of the zero-lift wing- 
fuselage pitching -moment coefficients from the A coefficient of equa- 
tion (B2) (for method I) and from the 6 q coefficient of equation (B5) 

(for method II) is illustrated with results presented for both methods. 

Method I .- Equation (A9) of the appendix is used to determine the 
wing -fuselage zero-lift pitching-moment coefficient as 


where 


C - Ax t 
Lm °M qSc 

*t = l t + d 


( 10 ) 


(ID 


For the example maneuver (flight 12, run 27 ) , the substitution into 
equation (10) of numerical values for A from equation (B4) and 

x t = (-552) + (-1.97) = -55^ inches 
produces a measured C~ value equal to 
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(-1T02) (-554) 

Cm 0M (159) (1428) (155.9) 


-0.0266 


The values of 


Cm 0M 


thus obtained from the 


A 


coefficients of table III 


and the use of equations (10) and (ll) are given in table VI and are 
plotted as a function of Mach number in figure 9 for each of the 68 maneu- 
vers. It was demonstrated in reference 3> however, that the measured tail 
loads were subject to large zero shifts (ranging from 16,000 to 
-6,000 pounds). In order to correct the measured zero-lift tail load, 
use Is made of the equation 


A s A - Z' - — 

A corr xt 


( 12 ) 


where Z' is the zero shift in measured tail load given in table III of 
reference 5 for the same maneuvers used in the present analysis and Zq> 
is the zero shift in tail aerodynamic torque determined from a tail torque 
and tail angle-of-attack analysis similar to that used for the tail loads 
in reference 5. For the example maneuver, use of equation (12) gives 


A cor r = ~ 1 T02 - (20) - (240) = -1962 

The zero-lift wing-fuselage pitching-moment coefficient is then recalcu- 
lated by use of the equation 




-A 


corr 




qSc 


(13) 


which with numerical values inserted becomes 


_ -(-1962) (-55*0 
Cm o " (159) (222625) 


-0.0307 


for the example maneuver. The error in may be estimated by the 

use of the standard error in the A coefficient as 


= £ A x t 
Cm 0 I qSc 

Using the standard error of A from table III gives 

ec = (*563) (55M_ = ±0 .0057 

“oi (159) (222625) 


(i*0 
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The corrected zero-lift wing-fuselage pitching-moment coefficients and 
their associated errors are listed in table VI. Figure 10 is a plot of 
the CmQ values as a function of Mach number and, although scatter does 

still exist, these results are a decided improvement over the results pre- 
sented in figure 9* 


Method II ,- The elevator-angle type of solution may be used to deter- 
mine the zero-lift wing-fuselage pitching-moment coefficients, as illus- 
trated by equation (A24) of the appendix: 


r 





3.00 + g — ^ 


t. de 
V 2 da 



( 15 ) 


The substitution of numerical values for the example maneuver, using the 
6 q value of table IV, geometric parameters, and tail angle -of -attack 
parameters from reference 3> results in 


^moji “ -0.0326 


with an associated error of ±0.0033. Values of and their standard 


errors are listed in table VI. Figure 11 is a plot of the C] 


mon 


value s 


as a function of Mach number, which is seen to indicate (although with 
somewhat more scatter) the same average variation with Mach number as 
illustrated in figure 10 for the method I data. 


The correlation between these two methods of evaluating zero-lift 
wing-fuselage pitching-moment data is seen more clearly in figure 12 
where Cm^ is plotted against The solid line in this figure 

is the perfect agreement line and the dashed sidebands represent an 
average departure from agreement based on the average of the errors 
listed in table VI for each method. With a few exceptions most of the 
data lie reasonably close to the correlation line. 

Variation with Mach number .- The CmQ data of figures 10 and 11 

indicated a tendency for the lower altitude data to have smaller absolute 
values of Cm q. This trend is in agreement with theory since increasing 
the dynamic pressure at constant Mach number relieves both the bending 
and torsional moments associated with the zero-lift wing loads and thus 
reduces the wing contribution to C m Q. Theoretical calculations indicate, 
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however, that for the maximum qm^ of the present tests the change in 
CmQ would be only 0.0010. Other factors not specifically corrected for 
in either set of data are the pitching-moment coefficient due to 

tail drag which was estimated to have a maximum value of 0.0023 and the 
pitching-moment coefficient due to engine thrust which was estimated to 
vary from 0.0019 to 0.0014 for the dynamic pressure range of the 68 maneu- 
vers used in the analysis. 

In order to determine a more definite variation of C m Q with Mach 

number, the data shown in figures 10 and 11 were used to determine weighted 
average values of at each of the Mach number groups shown in table VI. 

The weights were assigned by use of a weighting factor w similar to that 
given in equation (j) except that e X ac was replaced by ^he 

weighted average CmQ is given by the equation 



and the standard error of the weighted average is given by the equation 



The results of the application of equations (16) and (17) to the data 
for each of the Mach number groups are given in table VI and plotted in 
figure 13 as a function of Mach number. 

Comparison with wind-tunnel data .- Wind-tunnel data corrected for 
model flexibility effects are given in reference J on pages L-124, L-126, 
and L-128. From these data the wing-fuselage zero-lift pitching -moment 
coefficients (shown as the diamond symbols in fig. 13) were obtained. 


Pitching Moment of Inertia and Radius of Gyration 

In the following sections the effects of wing flexibility on the 
measured and calculated effective airplane moments of inertia are presented. 

ky^ from theory.- The effective moment of inertia of the airplane 
including the effects of wing flexibility is defined by equation (A10) as 
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d© 


( 18 ) 


The values of dLy/d© and dg to be used in equation (18) are difficult 

to determine by direct experimentation and recourse is made here to esti- 
mating them by theoretical means. The aerodynamic wing load due to 
pitching acceleration dL^/d© and the associated chordwise center of pres- 
sure dg were computed by the superposition method of reference 8 as modi- 
fied and used in references 4 and 5* The results of these computations are 
shown in figure 14 where the load (in this case for both wings) and the 
center of pressure (given in terms of percent M.A.C.) are plotted as func- 
tions of the flexibility parameter qm^. The relationship between the 

ordinate Xq of figure 14(a) and the center of pressure dgf is given 
by the equation 


d 8 - - X il)l5o 


(19) 


where X ac may be found from figure 8. The substitution of equation (19) 
into equation (18) and division by W/g produces the following equation 
for effective radius of gyration squared: 


k 2 - v 2 - § 
*Vf ■ *y w 


( x ac ~ x ©)^- 


0 /ioo d'e 


( 20 ) 


It was found that in the range of interest of qmp for the present tests 
the product ^X^ - Xg^^p was approximately, linear and was equal to 


(Xac - *e)j3* - 10,000qm R 


( 21 ) 


Thus, ky f in units of square feet is given by the equation 


k y 2 = k 2 - 32.2 x _ 1 .999 10 oo 0qm 

yi y 12W K 


tyf 


“ ky 2 - 1*2,000 2”E 


( 22 ) 


Values of ky f 2 calculated by equation (22) are given in table VII and 
are plotted in figure 15 as a function of qm^/W for the 68 test maneu- 
vers. In the calculation of ky 2 the empty-weight moment of inertia 
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given in reference 9 for a ground determination of Iy for the test 
airplane of 933,000 slug-feet 2 was used along with the fuel weights and 
locat ions . 


k v 2 from flight data.- The squared values of radius of gyration, 

y f* 

obtained by use of method I and method II flight data and equations (A10) 
and (A2 6 ) of the appendix (along with their standard errors), are given 

in table VII for all maneuvers. The average values of ky f 2 (from 


methods I and II) are also given in table VII and plotted in figure 15 • 
It will be noted that the measured k yf . 2 values have more scatter but 
a greater mean variation with qm^ than the calculated values. The 

Q 

disagreement between measured and calculated ky values may be due to 

actual differences in the rigid-wing values or an incorrect theoretical 
variation with qm^. In order to allow for these differences, the fol- 
lowing procedure was used to correlate flight and calculated values of 
ky^. 2 . If the values of ky^ 2 are assumed to be linear with respect to 

qmp, the following equations may be written: 


( k yf 2 )caic _ ( k y 2 + €l ) 


( dk yf 2 )calc q^R 

d q?°R w 

w 


(23) 



(24) 


Subtracting equation (23) from equation (24) results in the fol- 

2 

lowing equation for Aky^ : 
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Equation (25) now permits the correlation of all the data as a function 
of qniR since the actual ky 2 values that differ from run to run 

because of weight-distribution differences have been eliminated by the 

2 qm-n J 

subtraction. Columns of Ak v « values and — ~ are shown in table VII. 

A least-squares procedure was used with the data to establish the fol- 

2 

lowing relationship between Aky^ and qm^/W: 

Ak yf 2 = (6.6 ± 5.3) - (11.5 ± 2.6) X 10^ (26) 

Since the standard error of the value (€2 - €q) of equation (2 6) is ±5*5 
and the average standard error in the measured ky.^ from table VII is 
±5.8, the conclusion to be reached is that the flight measurements agree 
with the calculated values of ky 2 for rigid-wing conditions. The dif- 
ference between the measured and calculated k v 2 values is in the varia- 
qmo 

tion with -^=1. The equation which represents the flexible-wing effective 
radius of gyration is, thus, from equations (26), (25), and (22): 

ky f 2 = k y 2 - (115,000 ± 26,000 + 42,000)^p 

or 

ky f 2 = ky 2 - (157,000 ± 26,000)^ (27) 

2 

ky^ from wing -twist data.- The differences between the measured 

and calculated variations of ky^^ with qm^/W were sufficiently large 
to require some further evidence or confirmation. In figure 15 of refer- 
ence 2, optigraph measurements for the test airplane were reported which 
showed a considerable disagreement between measured and calculated values 
of wing twist associated with pitching accelerations. These twists as 
plotted in figure 15 of reference 2 are not those due to pitching inertia 
alone but include an additional air load component. After correction for 
this component, the aerodynamic load distribution due to twist resulting 
from pitching inertia was computed by use of the method of reference 10 
for the example of figure 15 in reference 2. Integration of the resulting 
load distribution for both wings gave a load per unit pitching accelera- 
tion of -15,040 pounds with a center of pressure at 86.9 percent c. 

Both of these points are indicated on the theoretical curves presented in 
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figure 14. Equation (20) was then used to compute the effective radius 
of gyration squared as 

k yf 2 = 358.35 - p86(1 - — — [(18.8 - 86.9) (-15,040)] 

** 144(108,000) L J 

= 318.52 feet 2 

where 358.35 was the rigid value of k y 2 and 108,000 was the weight for 
the conditions of the data of reference 2. The value of qm R for this 
example is 19*6; thus, the application of equation (27) results in 


k yf 


2 


358.35 - 1^0 00 _ 19<6 

108,000 

329.86 feet 2 


The theoretical value is 


determined from equation (22) as 350. 73 feet^. 


The following table indicates the correlation between the various 
methods of determining the effective radius of gyration squared for the 
test airplane under the average conditions used: 


Method 

V 2 ' 

ft 2 

e ky f > 

ft 2 

Theoretical fen. (PP)) 

350.73 

329.86 

318.52 


Analysis of all 68 flight maneuvers (eq. (27)) • • • 

Ui nor—'h'uri sh. mpfl filiT’pmpnt.fl 

±4.72 




0 

Better agreement is seen to exist between the ky f value using wing- 

twist measurements and the value determined by use of equation (27) than 
for the theoretical value determined by use of equation (22). 

DISCUSSION 


The preceding sections of the report have presented the results of 
the analyses and comparisons with theory and available wind-tunnel data 
as well as details of the analytical methods used to evaluate the flight 
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data. The following discussion concerns the importance of the many fac- 
tors included or omitted- in the present analysis. 

It is possible to lose sight of the importance of the supporting 
experimental studies necessary for the analysis of the horizontal-tail 
loads on a large flexible swept-wing aircraft. Although the present 
analysis was complicated by some difficulties, which were unforeseen 
when the flight-test program was originally laid out, the detailed sup- 
porting instrumentation which measured airplane angle of attack, wing 
deflection and twist, and fuselage deflection proved to be invaluable in 
many phases of the data analysis. In brief, the determination of the 
rigid wing-fuselage aerodynamic-center position from the flight data was 
considerably simplified as a result of the flight measurements of wing 
center of pressure presented in reference 5* 

The determination of wing-fuselage zero-lift pitching-moment coef- 
ficients from the flight data, which included large zero shifts, was only 
possible as a result of the tail-load analysis presented in reference 3 
which required the wing -lift -curve slope and angle-of-zero-lift data pre- 
sented in reference 4. A major factor in rationalizing the tail loads as 
functions of angle of attack and elevator angle was the availability of 
flight measurements of fuselage deflections which have since been reported 
in reference 11. 


The evaluation of the effects of wing flexibility on the effective 
moment of inertia was aided by the theoretical studies used in refer- 
ences 1, 2, 4, and the basic moment-of -inertia data provided by refer- 
ence 9, and the supporting check information obtained from the wing-twist 
measurements reported in reference 2. 

The analysis of elevator-angle -deflect ion time-history data used in 
the present report to confirm the direct tail-load evaluation of 
aerodynamic-center position, zero-lift wing-fuselage pitching-moment 
coefficient, and effective airplane moment of inertia was based directly 
on the analysis of references 3 and 4. 

Basic-Data Coefficients 


When method I was used to fit the tail-load time-history data, 
it was found that the wing-fuselage pitching moment due to pitching 
velocity produced immeasurably small horizontal tail loads. A theoretical 


02 t , 

calculation of the tail load per unit — ^4- (as used in the airplane 

pitching-moment equation (Al) , not the tail- load — angle-of -attack equa- 
tion (A15)) for the example maneuver (flight 12, run 27) resulted in 
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ez 


-18.0 pounds per degree 


dL tm 

which may be compared with the discarded value of — ^=— = 82 pounds 

0 - 1 + 


per degree noted for this maneuver in the flight analysis equation (B3). 


The equations used to fit the tail-load and elevator -angle time- 
history data (eqs. (B2) and (B5)) did not allow for the effects of dynamic 
wing and fuselage flexibility and, yet, no discernible differences were 
found in comparing the coefficients of the slow-rate and fast -rate maneu- 
vers, In several cases , dynamic wing-flexibility effects were suspected 
and the use of wing-tip flapping accelerations improved the fit to the 
time-history data without altering the primary coefficients. It was felt, 
therefore, that reasonably accurate results were thus obtained for the 
bulk of the data analyzed. 


A factor somewhat more difficult to define positively is the linearity 
of the wing-fuselage pitching-moment coefficient with normal-force coeffi- 
cient. The forms of the equations used force linearity and, for some of 
the higher altitude data where results were obtained at relatively high 
normal-force coefficients, this may have resulted in erroneous slopes and 
intercepts. All suspected departures from linearity were checked by the 
error-distribution time histories. 

Another probable source of error was the sloshing motion of the fuel 
in the three large unbaffled fuselage fuel tanks, which may have introduced 
errors in the assumed equation of airplane motion. 


Aerodynamic -Center Position 

Aside from the factors previously mentioned the accuracy with which 
the in-flight center-of-gravity position could be determined governs the 
accuracy of the flight values of wing-fuselage aerodynamic-center posi- 
tion. Corrections were made to account for the effect of airplane atti- 
tude on the fuel level and the resultant effect on center-of-gravity 
position. However, there was some indication from ground tests that the 
fuel-gage readings were not entirely independent of airplane attitude. 
The agreement shown between the aerodynamic -center positions determined 
by methods I and II is excellent and with few exceptions well within 
their calculated standard errors. The agreement between the data for 
the rigid-wing aerodynamic -center position shown in figure 7 and the 
wind-tunnel data seems reasonable. In the Mach number range from 0.70 
to 0.775 there is a difference of only 0.01c. The rearward shift of 
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aerodynamic-center position indicated above a Mach number of about 0.78 
is believed to be associated with the fuselage since the wing aerodynamic- 
center position given in reference 5 remains at roughly 23.4 percent of 
the local wing chord up to M = 0.8l. 


Zero-Lift Pitching -Moment Coefficient 

The agreement between the corrected zero-lift pitching-moment coef- 
ficients measured by use of method I and those measured by use of method II 
indicated that the tail-load zero-shift correction method proposed in ref- 
erence 3 was theoretically and practically sound. The differences between 
the flight values of shown in figure 13 and the wind-tunnel values 

(ref. 7) are large and not easily explainable. 

A source of error in the determination of Cjuq from flight data by 

use of either method I or method II is the accuracy of the elevator-angle 
measurement . While, for the present tests, the elevator-angle data were 
repeatable to within ±0.1 degree for comparable flight conditions, an 
error in CniQ due to elevator-angle zero errors /e c \ would be approxi- 
mately equal to \ ^B/ 


ep * 0.68Ct~€ r 
c ra 0& ^6 6 

which with a maximum value of Cj^ from reference 4 as 0.03 would be 



0.02e & 


As detailed in the section entitled "Methods and Results," the neg- 
lect of the effects of wing twist, tail drag, and pitching moment due to 
engine thrust could produce a maximum error in the flight measurements 
of CmQ of +0.0052. 


Moment of Inertia and Radius of Gyration 

The major factor affecting the accuracy of the determination of Iy f 
or for any individual maneuver is believed to be the effect of 

fuel sloshing on the airplane motion. In some abrupt maneuvers during 
which the tail of the airplane accelerated from positive to negative g 
and back again, the fuel in the rear tank slammed back to the bottom of 
the tank and, as a result, produced considerable vibration in accelerometer 
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readings at the center of gravity. Although an attempt was made to elimi- 
nate portions of the maneuvers from the analysis when such occurred, the 
large surging actions of the fuel undoubtedly have contributed to some 
of the scatter in the Iy data. 

The good agreement between the measured and calculated values of 
k v 2 for the rigid wing is believed to be indicative of the overall accu- 

racy of the methods used to extract the airplane static stability param- 
eters from the flight tail-load and elevator-angle measurements. 

The variation shown between theory and experiment in the changes 
of ky f 2 with the flexibility parameter qmp/W is puzzling even though 

it has been confirmed qualitatively by the wing-twist measurements of 
reference 2. This disagreement may be associated with unaccounted for 
local wing-section distortions near the tip or wing twists associated 
with dynamic motions of the wing tip which may be closely phased with 
the pitching acceleration time history. 


CONCLUDING REMARKS 

Flight measurements of wing-fuselage aerodynamic -center position, 
zero-lift pitching-moment coefficient, and effective airplane moment of 
inertia have been presented as derived from the analysis of 68 push-pull 
maneuvers on a large flexible swept-wing airplane in a Mach number range 
from 0.42 to 0.8l at pressure altitudes from 20,000 to 35>000 feet. The 
parameters as derived by two methods from measured horizontal -tail loads 
and elevator angles were in excellent agreement. The method for cor- 
recting for tail-load zero shifts (proposed in NACA RM L56J02) was applied 
to the flight data with good results, as evidenced by a comparison of the 
uncorrected and corrected zero— lift pitching— moment coefficients derived 
from the tail-load measurements. 

The effects of wing flexibility on the aerodynamic-center position 
and on the zero-lift pitching-moment coefficient were predictable by 
theory. The relief provided by the wing bending and twisting due to 
pitching— acceleration inertia loads to the effective moment of inertia 
was not predictable by theory. For the one check case available, essen- 
tial agreement was obtained between optigraph data and tail- load data 
for the relief due to pitching-acceleration wing loads. 

Specifically for the test airplane it was found that: 

1. Reasonable agreement was obtained between f light -measured rigid- 

ving fuselage aerodynamic-center positions and those determined from 

wind-tunnel tests. 
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2. The rigid-wing — fuselage aerodynamic -center position was con- 
stant at 24.4 percent mean aerodynamic chord up to a Mach number of O .78 
after which it shifted rearward as the Mach number increased to its maxi 
mum flight value of 0 . 8 l. 

5. The rearward shift in the wing-fuselage aerodynamic center above 
a Mach number of O .78 appeared to be associated with the fuselage since 
the wing aerodynamic-center position given in NACA RM L57E28 remained 
constant at 23.4 percent of the local wing chord up to a Mach number of 
0 . 81 . 


4. There were large differences between the wind-tunnel and flight - 
determined values of zero-lift wing-fuselage pitching-moment coefficient 
over the complete Mach number range of the tests. 

5* The moments of inertia determined from the flight measurements 
agreed with calculated values for the rigid-wing case. 

6 . Calculated effective moments of inertia, which included the 
relief due to wing bending due to pitching-acceleration inertia loads, 
did not agree with those measured by the present flight tests. 

7. For the test airplane, the pitching moment due to pitching- 
velocity loads on the wing was shown to be insignificant both theoreti- 
cally and from the tail-load measurements . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., August 1 , I 957 . 
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APPENDIX A 


PITCHING-MOMENT EQUATIONS 


The pitching moments on the airplane may he expressed by taking the 
summation of the moments about the wing-fuselage aerodynamic-center loca- 
tion as 


In equation (Al) positive forces act upwards, positive moments are nose 
up, and distances rearward of the aerodynamic center are negative. 

The parameters in equation (Al) which are functions of the wing 
flexibility are the aerodynamic-center location d, the wing-fuselage 
zero-lift pitching-moment coefficient C^, the lift on the wing due to 

the load induced by wing deflection due to pitching inertia Lq, and 
the lift on the wing due to the distribution of lift due to pitching 
velocity Lq. The term Lq is defined as 


I 


M ac = 0 = -nWd + C^qSS - I y 0 + Ig d- + L‘ d* + L^ + C^SqS^ 


(Al) 



(A2) 


and the term Lg is defined as 


dlv §l t 


(A3) 


V 



Method I - Direct Tail-Load Measurement 


Equation (Al) may be written as 


L t x t + C m „6qS-t-C-t = - c mo q - SS + nWd + 
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For the test data analyzed in the present report, the horizontal-tail 
aerodynamic load and the tail aerodynamic torque 


T t - c mg 6( l s t^t 


(A5) 


were both measured by means of strain gagesj thus, equation (A4) may be 
used as 

T t 1 T t “ C ido ( 1 s5 , nWd L dL w d .\ 0 dL w ® 2 t ( A A) 

+ n + ~ l y it sir * tv (A6) 


In a form simplified for least-squares analysis, equation (A6) may be 
written as 


01 + 

L^,p = A + Bn + C0 + D — ^ — 


(A7) 


From the coefficients A, B, C, and D of equation (A7) and the corre- 
sponding terms of equation (A6), the pitching-moment parameters of inter- 
est may be defined as aerodynamic-center location d 


d 


W - B 


zero-lift pitching-moment coefficient Cn^ 


Cmo 


Ax t 

qSc 


effective airplane moment of inertia Iy^ 


(A8) 


(A9) 



ae 


d e = Cx t 


and the pitching moment due to pitching velocity Lgd£ 


(A10) 


ez+ 

L e d e = D -ir 


(All) 
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Method II - Pitching -Moment Parameters From 


Elevator -Angle Measurements 


For the determination of pitching-moment parameters from elevator- 
angle measurements, a more complicated series of equations must be 
written which also include the effect of fuselage deflection on the 
horizontal-tail angle of attack. Equation (Al) is rewritten as 


Lt - "CmQ 


qSc 


nWd 



de x t 


C m 6 &( l s tctl 

V *t 1 



de 

v 


(A 12 ) 


Now 


Lt 


= C 




(A13) 


and can be written specifically for the test airplane by using 

equations (4) and (6) of reference 3: 


a t = -3.00 + a w 




Ilk de (n - d) 
V da g v 2 


7 de _ 
t da 



dOf 


di+ 


di t 


5 + nt + La- 
dd “ 


dn. 


dLt 


(A14) 


By the substitution of equation (Al4) into (A13) and factoring out the 
tail-load terms, the tail-load equation may be rewritten as 


Lt = 


V qSt 


1 - c V s * 


-5.00 ♦ <*£ - g) ♦ £ » - 


de 

3a 


^t d€/ -] \ r ^t 

g — — (n - 1) + — n t 

r2 da 


Y c. u.a dn^ 

Hereafter^ the expression — 




di+ 


(A15) 


is defined as the effective 


1 - c ^ qS * 

tail-lift-curve slope 
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With the use of the additional equations 


n t = n - 1.3520 

(Al6) 

: °^adj + ep C % c + 

(A17) 

CN - nW - Iyf ® 
UN Ac ~ qS Xj-Sq 

(A18) 


where equation (Al6) is the tail load factor expressed in terms of the 
center-of -gravity normal load factor and pitching acceleration and equa- 
tions (A17) and (Al8) are used in reference 3 to define the wing angle 
of attack, equations (A15) and (A12) can he combined to give 


5 = &n + 


86 


®^t . 86 n . 86 2f 

1 rr- + r— n + -tt y 

01+ v 8n 89 


(A19) 


where 


db 


i 1 + i)K) 


dlv 


d e 


f 62+ x t qS t 

d — - 






da t . P c t 

+ Cm 6 5+ 


(A20) 


S c 

_r __E_ + 


-mo 


&0 = 


S t*t 


V 3 - 00 + g ^ lS)( Cl °t)f ~ ( x ' S)( c V)f (°°»dj + *“) 






(A21) 


db 

3n 


Wd 


lS t x t 


dn^ ( Cl <xt)f ~ S v 2 da( Cl ^t)f ( L ~ da)( Cl ^t)f a^qS_ 


( Cl «*t)f d^ + Cm 6 


(A22) 
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and 


L yf 


d8 

be 


9S t Xt 


+ 1.352 


(V) 


f dnt 


"l 1 


de\ ( Cl <^) 

daj eupX^Sq 


fhf 


(°V)f ST + Cm *i H 


(A23) 


From equation (A21) the zero-lift pitching-moment coefficient C mo may 
be derived as 


•'m r 


s t*t 

Sc 


( c v) f lr + ^ a ■ (v) 




1 - 


iA“o adJ + 2 - 75 


(A24) 


The aerodynamic -center location d may be derived from equation (A22) 
as 


MV): 


d = 


rs t+ ir i+8 ^! i + ' 1 

On d6 dn t y2 da 



w 

da/ 

a^qS 


+ Q ?T 


(A25) 


W 


- C 


& dot , di t , g It de , ( ± de\ W 

chi d6 dn-t y2 da \ da/apqS 


- rv)t 

and the effective moment of inertia Iy f , from equation (A23) as 


Mf 


9St x t 


1 + 


(l-^ 

\ da/ 


MV), 


as 

de 


K) 


, r v 

f d6 m 5 x-t 


1.352 


(V) 


dl t 
f <int 


J Sa-c 


(A26) 
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APPENDIX B 

ANALYSIS OF BASIC DATA 


As shown in appendix A the wing-fuselage aerodynamic -center posi- 
tion, zero-lift wing-fuselage pitching-moment coefficient, and effec- 
tive moment of inertia may be determined from the coefficients of equa- 
tions (AT) and (A19) • The following sections give the least-squares 
procedures used for the determination of these coefficients and illus- 
trative examples for both methods I and II. 


Solutions by Use of Method I 

For analysis of flight tail-load measurements for each maneuver, 
equations of the form of equation (A7) may be solved for the coeffi- 
cients A, B, C, and D as follows: 



Ltrp 


(Bl) 


In equation (Bl) the individual rows of the rectangular matrix and the 
column matrix Lt T represent simultaneous measurements of the indi- 
cated parameters at each of various times in a given push-pull maneuver. 
Equation (Bl) in the following equivalent form was used for the deter- 
mination of the A, B, C, and D coefficients for each of the 
68 maneuvers listed in table II: 
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N 

I ” 

Z 8 

Z "‘i' 

V 


E* 

I” 2 

Z 5n 

Z 

B 

r 

— 





D 

^ J 


I 5 

Z ns 

Z 82 

Z ( 8, */ v ) 


I 5l t/ v 

£ "( S H/' r ) 

Z 5 K/ v ) 

Z KM 


— -L 

' 2N ' 

-< 

Z ^ 

> 


^ Ltj9 

1 

Z 

J 


(B2) 


The determination of the A, B, C, and D coefficients for each 
maneuver by the use of equation (B2) requires time-history measurements 
of the following quantities: 


n 

0 


«l t /V 

Ltij 


normal load factor at center of gravity 
pitching acceleration at center of gravity 
tail angle of attack due to pitching velocity 

aerodynamic tail load plus tail aerodynamic pitching moment 


Figures 2 and 3 illustrate in time-history form the data used in the 
analysis of a typical maneuver (flight 12, run 27) • The pitching accel- 
eration 0 has been corrected for instrument frequency-response 
characteristics . 

The use of equation (B2) with the time-history data of figures 2 
and 3 resulted in the following equation for tail load for this maneuver; 


Lt T = -1722 + 497n - 241500 + 820Z t /v (B 3 ) 

The coefficient of /V was small in comparison with its standard 

error and the fit to the time history was not improved by the inclusion 
of this pitching -velocity parameter. Analysis of the other maneuvers 
also indicated that the parameter OT-t/V did no " t; contribute significantly 

to the pitching-moment tail-load equation. Thus, for this and subsequent 
maneuvers the A, B, and C coefficients to be presented were obtained 


34 


NACA TN 4191 


by using equations of the form of equation (B2) with the 02^/v term 

omitted. With the pitching-velocity term omitted, the equation deter- 
mined as fitting the time history of 1^, for flight 12, run 27 became 

Lt T = “1702 + 392n - 240599 (B4) 

The calculated tail load L^ip (eq. (B4) ) is compared in figure 5 with 

the measured time history from which it was derived. The fit is reason- 
ably good and this sample maneuver is representative of the worst rather 
than the best correlation obtained. 

The A, B, and C coefficients obtained from the least-squares 
analysis of each of the 68 maneuvers are given in table III. Also listed 
in the table are the standard errors of fit and the standard errors of 
the individual coefficients. 


Solutions by Use of Method II 

The elevator-angle equation (eq. (A19)) may be written in a form for 
least-squares solution of time-history data as follows: 




N 

z» 

Z s 

£ eit/v 

-1 

r Z> ' 

d6 

bn 

< 

► = 

Z” 

Z” 2 

I s " 1 

) (Sl t /V)n 

< 

Z 6 ° 

> 

be 


I s 

Z ° s 

z s2 ; 

> (ez t /v)e 


Z 6S 

bb 


Z §it/v 1 

) n(5lt/V) | 

y e(ez t /v) ^ 

> ( 0 H/v) 2 


£ &(9i t /v) 

d<3l t /v 

V. J 




In figure 4 the elevator angle is shown in time-history form for 
the example maneuver (flight 12, run 27). When the elevator-angle data 
are used in addition to the data shown in figure 3, the following equa- 
tion is obtained by solution of equation (B5): 


= 4.548 + 3*534 


9 I t 


- 6.620n - 21.5638 


& 


V 


(b6) 


NACA TN 4191 


35 


In this case the tern Glt/V is an important parameter and is retained 

for all solutions. A time history of the elevator angles computed by use 
of equation (b 6) is shown in figure 4. The agreement is reasonably good 
but again the example maneuver is not as good a fit as was obtained with 
the majority of the maneuvers analyzed. 

The elevator -angle coefficients of equation (B5) for all 68 maneuvers 
are listed in table IV along with their standard errors and errors of fit. 
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TABLE I.- AIRPLANE CHARACTERISTICS AND DIMENSIONS 


Horizontal tail: 

Total area, sq ft 268 

Span, ft 33 

Root chord, ft 11.42 

Mean aerodynamic chord, ft 8.58 

Distance from horizontal-tail 25 percent M.A.C. to wing 

25 percent M.A.C., ft 46.52 

Incidence angle, deg -0.25 

Sweepback (25-percent-chord line), deg 32.9 

Aspect ratio 4.06 

Taper ratio 0.425 

Airfoil section BAC 100 

Strain-gage reference station (percent semispan) 5*3 

Wing: 

Total area, sq ft 1,428 

Span, ft Il 6 

Mean aerodynamic chord, in 155*9 

Aspect ratio 9*^-2 

Taper ratio 0.420 

Incidence angle, deg 2.75 

Sweepback (25-percent-chord line), deg 55 

Airfoil section BAC 145 
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TABLE III.- MEASURED TAIL-LOAD COEFFICIENTS OBTAINED BY USE OF EQUATION (B2) 


Flight 

Run 

A 

e A 

B 

£ B 

C 

e C 

Standard 

error 

2 

27 

-1,264 

±184 

279 

±180 

-30,725 

±980 

104 


28 

-4,207 

±288 

920 

±180 

-28,944 

±861 

122 


29 

-6,553 

±312 

776 

±121 

-29,533 

±811 

78 

3 

ii 

-2,143 

±207 

-1,260 

±149 

-29,823 

±849 

95 


12 

-1,627 

±99 

-1,446 

±84 

-31,656 

±611 

123 


13 

-799 

±93 

-1,420 

±76 

-30,567 

±722 

89 


14 

385 

±81 

-1,575 

±106 

-31,408 

±578 

69 

4 

19 

-6,451 

±191 

1,108 

±77 

-24,466 

±424 

69 


20 

-4,272 

±206 

691 

±153 

-24,964 

±500 

95 


21 

-1,181 

±185 

91 

±202 

-26,844 

±625 

94 

6 

11 

-3,782 

±123 

-1,558 

±101 

-26,492 

±626 

85 


12 

-4,343 

±214 

-1,016 

±92 

-22,263 

±485 

76 


13 

-3,352 

±357 

-1,305 

±211 

-23,929 

±755 

146 


14 

-2,062 

*545 

-1,610 

±229 

-26,189 

±984 

l4l 


15 

-1,494 

±242 

-1,286 

±182 

-24,484 

±793 

105 

8 

4 

8,375 

±292 

1,302 

±111 

-27,745 

±743 

69 


5 

6,254 

±132 

1,805 

±99 

-27,986 

±460 

101 


6 

2,096 

±50 

2,238 

±88 

-26,878 

±411 

89 

9 

1 

14,921 

±352 

1,259 

±144 

-31,453 

±553 

91 


2 

12,583 

±597 

1,553 

±121 

-30,431 

±1,012 

110 


3 

10,074 

±338 

1,526 

±114 

-30,379 

±670 

73 


4 

7,494 

±217 

1,948 

±144 

-28,856 

±493 

98 


5 

5,910 

±125 

2,073 

±79 

-27,788 

±46l 

89 


6 

2,723 

±190 

1,418 

±211 

-30,925 

±1,204 

108 


7 

1,815 

±187 

1,000 

±205 

-30,484 

±1,321 

91 

10 

3 

11,828 

±326 

988 

±218 

-30,098 

±727 

153 


4 

9,410 

±281 

1,091 

±285 

-32,209 

±875 

203 


5 

7,705 

±200 

1,124 

±262 

-29,850 

±686 

128 


6 

3,218 

±84 

1,775 

±136 

-28,107 

±4l6 

98 


7 

1,170 

±106 

1,554 

±168 

-29,071 

±556 

87 


8 

-673 

±296 

1,807 

±381 

-27,558 

±960 

182 


9 

-1,116 

±126 

1,379 

±130 

-28,669 

±487 

106 

11 

11 

362 

±386 

164 

±507 

-26,600 

±953 

295 


12 

-532 

±155 

417 

±183 

-26,567 

±553 

146 


13 

-1,650 

±104 

669 

±104 

-26,062 

±4o6 

89 


14 

-3,797 

±209 

1,057 

±156 

-24,022 

±521 

161 


15 

-4,566 

±263 

1,009 

±171 

-24,700 

±624 

210 


16 

-5,821 

±221 

1,257 

±142 

-23,845 

±4o6 

148 


17 

-7,315 

±267 

1,382 

±143 

-22,164 

±4l6 

124 


24 

-257 

±286 

501 

±216 

-26,232 

±1,263 

108 

12 

6 

6,965 

±216 

-716 

±320 

-26,798 

±598 

180 


7 

5,332 

±259 

-378 

±363 

-25,781 

±683 

284 


8 

4,988 

±473 

-1,873 

±658 

-27,161 

±943 

280 


9 

2,387 

±203 

-592 

±297 

-25,096 

±607 

279 


10 

-286 

±697 

-715 

±874 

-23,892 

±1,313 

480 


11 

-1,028 

±656 

-1,022 

±750 

-24,865 

±1,268 

366 


12 

-2,766 

±679 

-1,283 

±456 

-24,518 

±838 

169 


17 

-4,802 

±512 

-1,450 

±395 

-26,425 

±642 

222 


18 

-5,436 • 

±338 

-1,588 

±24^ 

-26,273 

±467 

152 


19 

-7,554 

±309 

-604 

±165 

-24,821 

±377 

191 


20 

-8,431 

±402 

-650 

±217 

-24,479 

±474 

186 


21 

-9,518 

±432 

-651 

±226 

-24,257 

±491 

155 


22 

-9,504 

±510 

-538 

±256 

-23,682 

±581 

211 


23 

-10,676 

±321 

-263 

±135 

-23,659 

±346 

130 


24 

-7,060 

±305 

1,338 

±184 

-22,970 

±440 

158 


25 

-5,765 

±259 

1,027 

±173 

-23,239 

±380 

163 


26 

-3,765 

±166 

842 

±128 

-24,437 

±35^ 

115 


27 

-1,702 

±365 

392 

±358 

-24,059 

±637 

267 


28 

-880 

±46l 

401 

±478 

-24,726 

±821 

289 

16 

1 

2,268 

±80 

-824 

±114 

-26,414 

±347 

90 


2 

2,367 

±138 

-870 

±158 

-25,959 

±454 

135 


3 

3,050 

±l4l 

-1,263 

±169 

-27,119 

±557 

117 


4 

3,202 

±134 

-1,467 

±161 

-27,302 

±496 

90 


5 

3,5 88 

±97 

-1,973 

±113 

-28,371 

±419 

71 


6 

2,789 

±69 

-1,677 

±80 

-27,929 

±242 

70 

17 

5 

-5,376 

±246 

-166 

±126 

-22,942 

±457 

135 


6 

-3,392 

±272 

-384 

±188 

-24,246 

±561 

221 


7 

-1,383 

±436 

-1,260 

±393 

-26,152 

±932 

225 
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TABLE IV.- MEASURED ELEVATOR-ANGLE COEFFICIENTS OBTAINED BY USE OF EQUATION (B5) 


Flight 

Run 

5 0 

E5 0 

86 

£ &>/d(ei t /v) 

dn 

€ &>/dn 

& 

€ 86/d0 

Standard 

error 

2 

27 

3.502 

±0.125 

4.648 

±0.163 

-6.514 

±0.145 

-27.369 

±0.244 

0.077 


28 

3.456 

±.145 

4 .606 

±.276 

-4.467 

±.150 

-19.595 

±.240 

.099 


29 

3.114 

±.078 

5.282 

±.233 

-3.357 

±.085 

-15.547 

±.183 

.048 

3 

11 

3.105 

±.157 

4.573 

±.244 

-6.123 

±.186 

-19.835 

±.270 

.068 


12 

3.408 

±.124 

3.500 

±.371 

-6.680 

± .144 

-22.237 

±.157 

.137 


13 

3.397 

±.127 

4.013 

±.289 

-7.672 

±.132 

-23.061 

±.324 

.111 


Ik 

3-895 

±.111 

3.212 

±.159 

-9.815 

±.146 

-51.507 

±.232 

.072 

4 

19 

4.130 

±.06l 

4.031 

±.173 

-3-995 

±.065 

-14 .042 

±.220 

.055 


20 

4.074 

±.151 

3.599 

±.190 

-5.763 

±.162 

-19.839 

±.200 

.106 


21 

4.055 

±.273 

3.251 

±.185 

-8.062 

±•339 

-27.583 

±.318 

.112 

6 

11 

3-219 

+ .067 

4.848 

±.225 

-4.303 

±.074 

-12.627 

±.121 

.047 


12 

3-225 

+ .098 

4.881 

±.216 

-4.279 

±.099 

-12.335 

±.132 

.060 


13 

3.447 

±.153 

4.101 

±.191 

-4.768 

±.150 

-13.523 

±.145 

.083 


14 

3.827 

±.190 

4.080 

±.243 

-5.876 

±.192 

-17.037 

±.212 

.089 


15 

3.484 

±.234 

3.862 

±.252 

-6.426 

±.254 

-18.269 

±.286 

.106 

8 

4 

5.553 

±.108 

3.321 

±.186 

-6.334 

±.119 

-22.399 

±.270 

.063 

5 

4.232 

±.063 

3.829 

±.135 

-4.597 

±.063 

-17.372 

±.110 

.070 


6 

3.945 

+ .035 

4.700 

±.102 

-3.163 

+ .037 

-12.170 

±.073 

.045 

9 

1 

2.641 

±.177 

3.936 

±.164 

-7.225 

±.224 

-30.945 

±.306 

.108 


2 

3.491 

±.249 

3.702 

±.204 

-6.594 

±.308 

-25.691 

±.390 

.082 


3 

3.564 

±.094 

4.362 

±.156 

-5.625 

±.111 

-21.500 

±.271 

.066 


4 

3.074 

±.192 

3-724 

±.284 

-4.570 

±.232 

-20.438 

±.528 

.088 


5 

2.960 

±.089 

4.387 

±.170 

-3.304 

±.076 

-17.589 

±.166 

.101 


6 

3.157 

±.147 

4.682 

±.309 

-3.398 

±.178 

-17.549 

±.324 

•080 


7 

3.335 

±.245 

5.067 

±.465 

-3.630 

±.3H 

-17.199 

±.535 

.109 

10 

3 

3.494 

±.142 

3.900 

±.104 

-6.650 

±.145 

-26.204 

±.l64 

.105 

4 

4.445 

±.189 

3.851 

±.168 

-5.994 

±.213 

-24.175 

±.220 

.143 


5 

4.088 

±.198 

4.333 

±.178 

-5.034 

±.222 

-20.248 

±.196 

.091 


6 

3.728 

±.098 

4.536 

±.151 

-4.357 

±.105 

-17.556 

±.117 

.081 


7 

3.885 

±.107 

4.506 

±.159 

-3.788 

±.122 

-15.881 

±.133 

.051 


8 

3.889 

±.170 

4.925 

±.249 

-3.447 

±.193 

-14.758 

*.175 

.084 


9 

3.328 

±.079 

5-353 

±.185 

-3.085 

±.080 

-13.770 

±.104 

.068 

11 

11 

4.312 

±.504 

3.609 

±.274 

-6.959 

±.617 

-25.417 

±.385 

.268 

12 

4.500 

±.175 

3.550 

±.150 

-6.032 

±.210 

-21.805 

±•195 

.112 


13 

4 .281 

±.113 

3.581 

±.142 

-5-244 . 

±.135 

-18.964 

±.163 

.072 


14 

4.411 

±.110 

3.883 

±.159 

-4.602 

±.122 

-16.518 

±.137 

.102 


15 

4.348 

±.095 

4.251 

±.l6l 

-4.067 

±.106 

-15.095 

±.121 

.107 


16 

4.359 

±.119 

4.128 

±.180 

-3.777 

±.137 

-13.936 

±.128 

.106 


17 

4.162 

±.088 

4.199 

±.147 

-3.525 

±.105 

-12.152 

±.103 

.082 


24 

3.805 

±.357 

3.2&7 

±.360 

-7.096 

±.321 

-28.281 

±.416 

.103 

12 

6 

1.912 

±.488 

4.254 

±.256 

-8.187 

*•579 

-26.922 

±.339 

.210 

7 

2.396 

±.396 

4.254 

±•237 

-7.097 

±.448 

-23.442 

±.249 

.226 


8 

2.370 

±.428 

4.902 

±.197 

-5.917 

*•507 

-20.762 

±.278 

.156 


9 

4.031 

±.201 

4.124 

±.191 

-6.808 

*.259 

-19.980 

±.154 

.209 


10 

3.058 

±.540 

5.052 

±.468 

-4.926 

±.670 

-15.822 

±.335 

.346 


11 

3.744 

±.493 

4.701 

±•541 

-5.403 

*•594 

-16.056 

±.327 

.256 


12 

3.726 

±.324 

5.008 

±.272 

-4.991 

±.404 

-14.843 

±.239 

.132 


17 

3.622 

±.631 

3.611 

±.219 

-9.588 

±.694 

-28.126 

±.316 

.283 


18 

4.059 

±.328 

3.582 

±.165 

-8.592 

*•375 

-24.189 

*•199 

.151 


19 

4.503 

±.144 

3.663 

±.102 

-7.185 

±.134 

-20.376 

±.090 

.141 


20 

4.594 

±.141 

3.856 

±.094 

-6.500 

±.149 

-17.992 

±.095 

.114 


21 

4.632 

±.167 

4.053 

±.143 

-5.772 

±.194 

-15.667 

±.123 

.097 


22 

4.212 

±.186 

4.265 

±.178 

-5.282 

±.214 

-14.488 

±.142 

.137 


23 

4.252 

±.103 

4.915 

±.112 

-4.686 

±.107 

-13.443 

±.082 

.094 


24 

4.438 

±.094 

4.038 

±.104 

-3.892 

±.100 

-13.193 

±.082 

.091 


25 

4.378 

±.109 

3.936 

±.096 

-4.651 

±.110 

-15.534 

*.079 

.104 


26 

4.072 

±.104 

3.683 

±.100 

-5.202 

±.110 

-17.859 

±.098 

.093 


27 

4.548 

±.239 

3.534 

±.130 

-6.620 

±.258 

-21.563 

±.144 

.180 


28 

3.985 

±.468 

3.393 

±.191 

-7.772 

±.500 

-26.485 

±.267 

.278 

16 

1 

4.573 

±.078 

3.805 

±.087 

-5.475 

±.08l 

-15.842 

±.084 

.062 

2 

3-955 

±.111 

3.913 

±.106 

-5.883 

±.114 

-16.838 

±.092 

.081 


3 

4.196 

±•159 

3.561 

±.121 

-7.165 

±.170 

-20.882 

±.149 

.083 


4 

4.164 

±.271 

3.390 

±.163 

-8.688 

±.293 

-25.164 

±.241 

.104 


5 

4.023 

±.151 

3.388 

±.087 

-10.629 

±.165 

-30.696 

±.158 

.072 


6 

3.181 

±.159 

3.659 

±.085 

-9.254 

±.173 

-28.111 

±.137 

.107 

17 

5 

3.937 

±.089 

5.758 

±.155 

-3.569 

±.094 

-10.948 

±.100 

.086 

6 

3-777 

±.139 

4.699 

±.183 

-4.075 

±.144 

-12.114 

±.125 

.142 


7 

4.094 

±.261 

5.042 

±.231 

-4.721 

±.292 

-13.946 

±.191 

.115 
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TABI£ V.- AERODmMIC -CENTER POSITION 


Group 

M 

Group 

( X ac)r 

Flight 

Run 

M 

qm R 

( X ac)flex 
(method I) 

( X ac)flex 
(method II) 


( X ec)r t 
(method I) 

( X ac)r 
(method II ) 

0.429 

23.71 ± °-52 

11 

12 

16 

16 

24 

28 

1 

0.427 
.427 
.428 
• 433 

11.4 

11.4 

11.5 
11.9 

22.1 ± O.7 
21.8 ± 1.5 
21.3 ± 0.3 

20.1 ± 0.3 

20.6 ± 0.9 

20.7 ± 1.4 

19.8 ± 0.4 
16.4 ± 0.5 

3-5 

3.5 
3-5 

3.6 

25.6 

25.3 

24.8 

23.7 

24.1 

24.2 
23-3 
20.0 

.486 

24.45 ± 0.35 

12 

16 

12 

4 

11 

27 

4 

17 

21 

11 

.482 

.482 

.483 

.486 

•495 

14-7 

14.8 
12.0 

11.9 
12.8 

21.6 ± 1.2 

19.4 ± 0.5 

19.7 ± 1.2 

21.8 ± 0.7 

22.5 ± 1.6 

20.8 ± 0.9 
18.7 ± 0.9 
18.0 ± 1.9 

21.9 ± 1.0 

20.3 ± 1*9 

4.4 

4.4 

3.6 

3.6 

3-9 

26.0 

23-8 

23.3 
25-4 

26.4 

25.2 

23.1 
21.6 
25-5 

24.2 

.541 

24.10 ± 0.27 

12 

11 

16 

12 

8 

18 

12 

3 

26 

4 

• 532 
.542 
.542 
.543 
.544 

14.8 

15.6 

19.0 

19.2 

15.5 

19.7 ± 0.8 

21.1 ± 0.6 

18.1 ± 0.5 
19.5 ± 0.4 

L 9.9 ± 0.3 

18.8 ± 1.2 

19.7 ± 0.7 

16.8 ± 0.6 
18.4 ± 0.4 
L8.5 ± 0.4 

4.4 
4.6 
5-5 

5.5 

4.6 

24.1 

25-7 

25-6 

25.0 

24-5 

23-2 

24.3 

22.3 
23-9 
23-1 

• 595 

24.33 ± 0.31 

12 

4 

12 

11 

10 

9 

16 

12 

6 

20 

25 

13 

3 

1 

2 

19 

.384 

• 591 

• 595 

• 597 
.598 
.598 

• 599 
.600 

12.3 

18.5 

23.6 

18.9 

15.5 

12.2 

24.0 

19.3 

18.0 ± 1.0 
19.2 ± 0.4 

18.5 ± 0.6 

20.1 ± 0.3 

20.9 ± 0.6 

20.6 ± 0.4 
17.4 ± 0.5 

16.2 ± 0.5 

14.8 ± 1.6 

20.6 ± 0.6 

18.2 ± 0.5 

19.2 ± 0.5 

19.8 ± 0.4 

18.8 ± 0.6 

15.6 ± 0.5 

18.3 ± 0.5 

3-7 

5-4 

6.7 

5-5 

4.6 

U 

5.6 

21.7 

24.6 
25-2 

25.6 
25.5 
24.3 
24.2 

21.8 

I8.5 

26.0 

24.9 
24-7 

24.4 

22.5 
22.4 

23.9 

.635 

24.23 ± 0.33 

3 

2 

11 

12 

14 

27 

14 

20 

.631 

.636 

.637 

.637 

13.9 

13.6 

22.1 

22.2 

19.2 ± 0.3 

21.3 ± 0.6 
18.7 ± 0.5 
16.3 ± 0.7 

21.3 ± 0.4 

19.5 ± 0.4 

18.5 ± 0.5 
17.9 ± 0.6 

4.2 
4.1 

6.3 
6.3 

23.4 

25-4 

25.0 

22.6 

25*5 

23-6 

24.8 

24.2 

.644 

23.96 ± 0.21 

12 

12 

16 

6 

9 

10 

8 

£ 

1 

15 

2 

4 

5 

.642 

.642 

.642 

.643 

.647 

.647 

.648 

14.7 

18.7 
14-7 
18.5 

23-4 

17.0 ±1.1 
17.4 ± 0.6 

17.3 ± 0.4 
17.9 ± 0.6 

19.3 ± 0.3 

20.0 ± 0.8 

18.1 ± 0.3 

14.9 ± 1-4 
18.2 ± 0.5 
15.8 ± 0.4 
17.1 ± 1.0 
19.7 ± 0.9 
19-7 ± 0.7 
I7.O ± 0.2 

4.4 
7-7 
7*7 

5.4 

4.4 

5.4 
6.6 

21.4 

25.1 

25-0 

23.3 

23 - 7 

25.4 

24 - 7 

19-3 

23.9 

23-5 

22.5 

24.1 

25.1 
23-6 

.681 

24.25 ± 0.36 

12 

9 

10 

11 

12 

8 

3 

5 

15 

21 

.681 

.681 

.682 

US.6 

16.5 

20.5 
25.3 
26.2 

21.1 ± 1-9 

19.4 ±0.3 
19.9 ± 0.7 

19.1 ± 0.6 

16.5 ± 0.7 

13.0 ± 1.7 
18.4 ± 0.3 
17.8 ± 0.8 
18.2 ± 0.5 
17.7 ± 0.9 

4.9 

4.9 

5.9 
7-1 
7.2 

26.0 

24.3 

25-8 

26.2 

23-7 

17.9 

23.3 
23.7 

25.3 

24.9 

.695 

24.32 ± 0.28 

3 
6 

12 

4 
11 

13 

14 
22 

16 

.689 

.690 

.694 

.699 

.702 

17.2 

22.2 
27.3 
27.7 
28.0 

18.5 ± 0.2 
18.9 ± 0.8 

l£.0 ± 0.8 

17.5 ± 0.3 
17.7 ± 0.5 

18.6 ± 0.4 

17.7 ± 0.8 

15.9 ± 1.0 

17.9 ± 0.3 
17.0 ± 0.7 

5.0 

6.3 

7 *i 

7.6 

7-7 

23-5 

25.2 

23-5 

25.1 

25-4 

23.6 

24.0 

23.4 

25.5 
24-7 

.726 

24.44 ± O.17 

12 

17 

10 

3 

9 

9 

l 

12 

4 

.721 

.725 

.726 

.728 

•731 

19.3 
32.2 
25.1 
20.6 

20.4 

17.0 ± 0.9 

18.0 ± 1.2 
17.9 ± 0.4 
18.5 ± 0.3 
18.2 ± 0.4 

19.0 ± 0.9 
14.9 ± 1.4 
17.8 ± 0.4 
18.4 ± 0.5 

18.0 ± 0.8 

5.6 

8.5 

7.0 

5-9 

5.9 

22.6 

26.5 

24.9 

24.4 

24.1 

24.6 

23.4 

24.8 
24.3 

23.9 

• 735 

24.89 ± 0.35 

11 

12 

2 

6 

17 

23 

28 

13 

• 734 

• 735 
•735 
.741 

32.3 
33.1 

20.4 

27.4 

17.3 ± 0.5 
15.2 ± 0.4 
18.9 ± 0.6 
17.6 ± 0.7 

18.1 ± 0.5 
15.9 ± 0.5 
18.4 ± 0.6 

16.2 ± 0.7 

8.6 

8.7 

5.9 

7.5 

25.9 

23-9 

24.8 

25.1 

26.7 

24.6 
24.3 

23.7 

.758 

25-53 ± 0.25 

3 

8 

17 

10 

11 

6 

6 

7 

.750 

.758 

.762 

.763 

22.3 

36.3 
38.0 
29.6 

17.9 ± 0.4 

16.8 ± 0.3 
15.4 ± 0.6 

18.8 ± 0.5 

18.2 ± 0.7 
16.0 ± 0.2 
l 4.6 ± 0.7 
18.5 ± 0.5 

6.4 

9-4 

9.7 

8.0 

24.3 

26.2 

25.1 

26.8 

24.6 

25.4 
24.3 

26.5 

.776 

24.41 ± 0.14 

12 

9 

10 

5 

•773 

•779 

24.0 

25.7 

17.4 ± 2.6 

17.5 ± 0.2 

17.0 ± 2.8 
16.9 ± 0.3 

6.8 

7-1 

24.2 

24.6 

23.8 

24.0 

.791 

26.58 ± 0.20 

6 

10 

6 

12 

9 

2 

11 
8 

12 
11 

6 

29 

.789 
.789 
.79 0 
.790 
.795 
.796 

32.2 

32.8 

26.3 

26*7 

26.7 

18.4 ± 0.3 

18.2 ± 1.1 

16.7 ± 0.3 
18.0 ± 2.8 

19.7 ± 0.6 

19.2 ± 0.4 

18.4 ± 0.4 

19.2 ± 0.9 
18.6 ± 0.5 

21.2 ± 2.6 
19.0 ± 0.7 

19.2 ± 0.4 

8.6 

8.5 

8.7 

7 ‘? 

7 .4 

7.4 

27.0 
26.7 

25-4 

25.3 

27.1 

26.6 

27.0 

27.7 

27-3 

28.5 
26.4 

26.6 

.810 

27.76 ± 0.48 

17 

9 

10 

12 

5 

7 

9 

12 

.808 

.810 

.812 

.812 

46.0 

28.5 

34.9 

29.0 

14.6 ± 0.4 
21.1 ± 0.6 

19.6 ± 0.4 

18.6 ± 1.4 

15.7 ± 0.5 
22.0 ± 1.3 

19.7 ± 0.4 
21.4 ± 1.8 

ll.l 

7.8 

9.1 

7-9 

25*7 

28.9 

28.7 

26.5 

26.8 

29.8 
28.8 
29.3 
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TABLE VI.- ZERO-LIFT WING-FUSELAGE PITCHING-MOMENT COEFFICIENTS 


Group 

M 

Group 

Cm 0 

Flight 

Run 

M 

C m 0 

(method I) 

Cffl 0 

(method I 
corrected 
for zero shift) 

e Cmo 

(method I) 

C m 

®0 

(method II) 

€ Cmo 

(method II) 



11 

24 

0.427 

-0.0051 

-0.0285 

±0.0057 

-0.0402 

±0.0063 

0.429 

-0.0371 ± 0.0023 

12 

28 

.427 

-0173 

-.0313 

±.0091 

-.0379 

+.0083 



16 

5 

.428 

.0704 

-.0319 

±.0019 

-.0389 

±.0027 



16 

6 

.433 

• 0532 

-.0358 

±.0013 

-.05^6 

+.0028 



12 

27 

.482 

-.0266 

-.0307 

±.0057 

-.0326 

±.00 33 

.486 

-0.0363 ± 0.0012 

16 

4 

.482 

.0501 

-.0361 

±.0021 

-.0401 

±.oo46 



12 

17 

.483 

-.0924 

-.0387 

±.0098 

-.0448 

±.0114 



4 

21 

.486 

-.0227 

-.0415 

±.0035 

-.0397 

±.0049 



11 

11 

.495 

.0065 

-.0264 

±.0069 

-.0354 

±.0089 



12 

18 

.532 

-.0866 

-.0388 

±.0054 

-.0413 

±.0057 

.541 

-0.0426 ± 0.0019 

11 

12 

.542 

-.0081 

-.0301 

±.0024 

-.0358 

±.0030 



16 

3 

.542 

.0383 

-.0394 

±.0018 

-.0432 

±.0025 



12 

26 

.543 

-.0466 

-.0441 

±.0021 

-.0452 

±.0016 



8 

4 

.544 

.1286 

-.0449 

±.0045 

-.0515 

±.0019 



12 

6 

.584 

.1379 

-.0554 

±.0043 

-.0823 

±.0092 

•595 

-0.0459 ± 0.0016 

4 

20 

.591 

-.0563 

-.0532 

±.0027 

-.0454 

±.0025 



12 

25 

.595 

-.0597 

-.0452 

±.0027 

-.0435 

±.0016 



11 

13 

.597 

-.0212 

-.0415 

±.0013 

-.0429 

±.0018 



10 

3 

•598 

.1851 

-.0492 

±.00 51 

-.0534 

±.0025 



9 

1 

.598 

.2981 

-.0507 

±.0070 

-.0682 

±.0034 



16 

2 

.599 

.024 2 

-.0455 

±.0014 

-.0496 

±.0029 



12 

19 

.600 

-.0963 

-.0546 

±.0039 

-.0379 

±.0023 



3 

14 

.631 

.0069 

-.0595 

±.0014 

-.0465 

±.0020 

.635 

-0.0495 ± 0.0026 

2 

27 

.636 

-.0229 

-.0478 

±.0033 

-.0551 

±.0023 



11 

14 

.636 

-.0429 

-.0446 

±.0024 

-.0431- 

±.0017 



12 

20 

.637 

-.0954 

-.0485 

±.0045 

-.0392 

+.0022 



12 

7 

.642 

.0914 

-.0678 

±.0044 

-.0753 

+.0074 

.644 

-0.0445 ± 0.0011 

12 

24 

.642 

-.0638 

-.0457 

±.0028 

-.0448 

±.0013 



16 

1 

.642 

.0203 

-.0417 

±.0007 

-.0445 

±.0011 



6 

15 

.643 

-.0201 

-.0489 

±.0033 

-.0596 

±.0038 



9 

2 

.647 

.2108 

-.0597 

±.0102 

-.0562 

±.0044 



10 

4 

.647 

.1270 

-.0417 

±.0038 

-.0403 

±.0032 



8 

5 

.648 

.0674 

-.0436 

±.0014 

-.0472 

±.0010 



12 

8 

.679 

.0766 

-.0447 

±.0073 

-.0771 

±.0079 

.681 

-0.0478 ± 0.0018 

9 

3 

.681 

.1566 

-.0539 

±.0053 

-.0553 

±.0017 



10 

3 

.681 

.0962 

-.0419 

±.0025 

-.0478 

±.0033 



11 

15 

.681 

-.0462 

-.0441 

±.0027 

-.0462 

±.0015 



12 

21 

.68 2 

-.0942 

-.0470 

±.0043 

-.0418 

±.0025 



3 

13 

.689 

-.0120 

-.0614 

±.0014 

-.0585 

±.0023 

.695 

-0.0523 ± 0.0015 

6 

14 

.690 

-.0240 

-.0517 

±.0040 

-.0529 

±.0031 



12 

22 

.694 

-.0916 

-.0505 

±.0049 

-.0487 

±.0028 



4 

19 

.699 

-.0615 

-.0530 

±.0018 

-.0496 

±.0009 



11 

16 

.702 

-.0550 

-.0485 

±.0021 

-.0483 

±.0017 



12 

9 

.721 

.0338 

-.0615 

±.0029 

-.0514 

±.0036 

.726 

-0.0620 ± 0.0013 

17 

7 

.725 

-.0118 

-.0513 

±.0037 

-.0556 

±.0037 



10 

6 

.726 

.0352 

-.0612 

±.0009 

-.0593 

±.0016 



3 

12 

.728 

-.0217 

-.0651 

±.0013 

-.0633 

±.0022 



9 

4 

.731 

.1017 

-.0728 

±.0029 

-.0711 

±.0035 



11 

17 

.734 

-.0632 

-.0592 

±.0023 

-.0557 

±.0013 

.756 

-0.0586 ± 0.0016 

12 

23 

.735 

-.0908 

-.0598 

±.0027 

-.0547 

±.0015 



2 

28 

.735 

-.0573 

-.0632 

±.0039 

-.0645 

±.0026 



6 

13 

.741 

-.0345 

-.0654 

±.0037 

-.0674 

±.0024 



3 

11 

•750 

-.0275 

-.0777 

±.0027 

-.0725 

±.0028 

.758 

-0.0637 ± 0.0007 

8 

6 

.758 

.0168 

-.0634 

±.0004 

-.0628 

±.0005 



17 

6 

.762 

-.0263 

-.0664 

±.0021 

-.0653 

±.0019 



10 

7 

.763 

.0115 

-.0652 

±.0010 

-.0629 

±.0017 



12 

10 

.773 

-.0036 

-.0789 

±.0088 

-.0774 

±.0098 

.776 

-0.0809 ± 0.0003 

9 

5 

.779 

.0460 

-.0807 

±.0015 

-.0812 

±.0016 



6 

11 

.789 

-.0360 

-.0733 

±.0012 

-.0772 

±.0011 

•791 

-0.0770 ± 0.0010 

10 

8 

.789 

-.0065 

-.0733 

±.0029 

-.0659 

±.0028 



6 

12 

.790 

-.0428 

-.0791 

±.0020 

-.0772 

±.0016 



12 

11 

.790 

.0024 

-.0802 

±.0104 

-.0673 

±.0088 



9 

6 

• 795 

.0315 

-.0801 

±.0022 

-.0798 

±.0027 



2 

29 

.796 

-.0760 

-.0838 

±.0036 

-.0813 

±.0014 



17 

5 

.808 

-.0374 

-.O76O 

±.0017 

-.0700 

±.0012 

.810 

-0.0779 ± 0.0023 

9 

7 

.810 

.0201 

-.0866 

±.0021 

-.0785 

+ .0044 



10 

9 

.812 

-.0102 

-.0826 

±.0012 

-.0795 

±.0013 



12 

12 

.812 

-.0504 

-.0848 

±.0075 

-.0715 

±.0058 
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TABLE VII.- RAD IUS-OF -GYRATION DATA 


Flight 

Run 

lc 2 

k y f 

(method II) 

V 

(method II) 

(method I) 

V 

(method I) 

( k y/)av 

(methods I and II) 

(^f )calc 

l m R 

W 


2 

27 

28 
29 

370.8 

357.8 
348.5 

±3.2 

±4.0 

±3.7 

404.9 

385.9 
394.1 

±12.9 

±11.5 

±10.8 

387.8 

371.8 
371.3 

363.0 

360.7 

358.3 

1.208 X 10- 4 
1.834 
2.380 

24.8 

ll.l 

13.0 

3 

11 

12 

13 

14 

359.1 
377.3 
351.6 

404.2 

±4.5 

±2.5 

±4.6 

±2.8 

371.8 

394.7 

382.8 
395.2 

±10.6 

±7.6 

±9.0 

±7.3 

365.4 

386.0 

367.2 

399.7 

367.5 

367.9 
369.2 

369.9 

1.854 

1.715 

1.435 

1.168 

-2.1 

18.1 

-2.0 

29.8 

4 

19 

20 
21 

317.0 

341.5 

340.8 

±4.5 

±3.2 

±3.8 

337.6 

343.9 

367.5 

±5.9 

±6.9 

±8.5 

327.3 

342.7 

354.2 

335.1 

338.3 

340.6 

2.544 

1.702 

1.098 

-7.8 

4.4 

13.6 

6 

11 

12 

13 

14 

15 

318.0 

314.1 

306.2 
340.3 
313.5 

±2.7 

±3.0 

±3.0 

±3.9 

±4.6 

365.2 
308.9 
331.7 
362.4 

341.3 

±8.6 

±6.8 

±10.4 

±13.6 

±11.0 

341.6 

311.5 

319.0 

351.3 

327.4 

343.7 

343.8 
345.6 
348.1 
349.3 

2.960 

3.017 

2.528 

2.052 

1.738 

-2.1 

-32.3 

-26.6 

3.2 

-21.9 

8 

4 

5 

6 

304.3 

315.3 
290.7 

±3.4 

±1.8 

±1.5 

332.3 
337-2 

326.3 

±8.9 

±5.6 

±5.0 

318.3 

326.2 

308.5 

351.7 

349.4 

345.4 

1.242 

1.880 

2.927 

-33.4 

-23.2 

-36.9 

9 

1 

2 

3 

4 

5 

6 
7 

344.3 

335.0 
307.7 

337.1 
344.6 

348.5 

360.6 

±3.2 

±4.7 

±3.6 

±4.9 

±2.9 

±5.8 

±9.9 

370.4 

360.4 
360.1 
344.3 
333.0 

368.5 
362.8 

±6.2 

±12.0 

±8.0 

±5.9 

±5.5 

±14.4 

±15.7 

357.3 

347.7 
333.9 

340.7 

338.8 
358.5 
361.7 

363.9 

362.7 
362.1 
361.2 

359.8 
359.7 
359.1 

.963 

1.165 

1.308 

1.623 

2.049 

2.133 

2.282 

-6.6 

-15.0 

-28.2 

-20.5 

-21.0 

-1.2 

2.6 

10 

3 

4 

5 

6 

7 

8 
9 

347.0 

564.0 
334.8 
329.6 

332.3 
328.2 

327.4 

±2.0 

±3.1 

±3.0 

±2.0 

±2.4 

±3.4 

±2.1 

351.8 

380.0 
352.6 

334.9 

347.1 

330.2 

343.1 

±8.5 

±10.3 

±8.1 

±4.9 

±6.4 

±11.5 

±5.8 

349.4 

372.0 

343.7 

332.2 

339.7 

329.2 

335.2 

360.5 

359.7 
359.3 
357-8 
356.9 

356.5 

355.8 

1.219 

1.462 

1.623 

1.990 

2.360 

2.532 

2.794 

- 11.1 

12.3 

-15.6 

-25.6 

-17.2 

-27.3 

-20.6 

11 

n 

12 

13 

14 

15 

16 
17 
24 

336.1 
534.6 

536.5 

327.2 

328.6 

519.2 
503.7 

355.2 

±4.8 

±2.8 

±2.7 

±4.3 

±2.4 

±2.6 

±2.3 

±5.0 

361.5 

362.7 

358.4 

331.5 

343.8 

332.4 

309.8 

375.4 

±13.0 

±6.6 

±5.6 

±7.2 

±8.7 

±5.7 

±5.8 

±18.1 

348.8 

348.6 
347.4 
329.3 

336.2 

325.8 

306.7 

365.3 

356.9 

356.0 

355.7 

354.5 
353.3 

352.5 

350.5 

359.5 

1.172 

1.433 

1.742 

2.037 

2.334 

2.597 

3.005 

1.099 

-8.1 

-7.4 

-8.3 

-25.2 

-17.1 

-26.7 

-43.8 

5.8 

12 

6 

7 

8 
9 

10 

11 
12 

17 

18 

19 

20 

21 
22 

23 

24 

25 

26 

27 

28 

318.0 

519.1 

314.9 

326.9 

308.6 

330.1 

327.9 

531.6 

334.1 
340.5 

530.4 

521.3 

307.7 

313.5 

294.3 

309.9 
312.0 

311.1 

314.1 

±3.8 

±3.2 

±2.5 

±2.3 

±6.0 

±6.1 

±4.8 

±3.6 

±2.6 

±1.4 

±1.1 

±2.3 

±2.7 

±1.7 

±1.6 

±1.5 

±1.6 

±2.0 

±3.0 

334.4 

329.7 

336.7 
315.9 
302.0 

314.0 

309.3 

338.7 

337.1 

323.3 

317.4 

315.1 

308.4 

310.5 

311.2 

313.2 

330.2 

324.2 
333.0 

±7.5 

±8.5 

±11.9 

±7.7 

±16.6 

±7.3 

±10.6 

±8.2 

±6.0 

±5.2 

±6.1 

±6.4 

±7.6 

±4.5 

±6.0 

±5.1 

±4.8 

±8.6 

±11.0 

326.2 

324.4 

325.8 

321.4 
305.3 
322.0 
318.6 

335.2 

335.6 

331.9 

323.9 

318.2 
308.0 
312.0 

302.7 

311.5 
321.1 

317.6 
323.5 

350.1 

349.3 

348.7 

347.4 
346.3 

345.5 

344.2 

349.3 

548.4 

347.1 

546.1 

344.7 

544.8 

345.4 

336.2 

337.8 

339.5 

341.3 

541.8 

1.022 

1.222 

1.384 

1.614 

2.015 

2.214 

2.443 

1.029 

1.270 

1.658 

1.909 

2.257 

2.358 

2.868 

2.511 

2.124 

1.736 

1.333 

1.034 

-23.9 

-24.9 

-22.9 

-26.0 

-41.0 

-23.5 

-25.6 

-14.1 

-12.8 

-15.2 
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Figure 1 .- Side view of test airplane. L- 86692 
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Figure 2.- Time-history data for push-pull maneuver. 
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Figure 3.- Time histories of measured and calculated horizontal- tail aerodynamic loads for 

maneuver of figure 2. 
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Figure 4.- Time histories of measured and calculated elevator angles for maneuver of figure 2. 
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Figure 5.- Forward shift in wing aerodynamic-center position as a func- 
tion of flexibility parameter from figure k(c) of reference 5* 
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Figure 7.- Comparison of flight and wind-tunnel wing-fuselage aerodynamic-center positions for 

rigid-wing, conditions. 
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Figure 8.- Wing-fuselage aerodynamic -center position for flexible-wing conditions. 
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Figure 9.- Zero-lift wing-fuselage pitching-moment coefficient measured 
by method I as a function of Mach number. 
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Figure 10.- Zero-lift ving-fuselage pitching-moment coefficients as obtained by method I and 

corrected for tail- load zero shifts. 
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Figure 11.- Zero-lift wing-fuselage pitching-moment coefficients as obtained by method II. 


¥ 


I 

> 


■r 

H 

VO 


6tj AJef-/io<y T 


NACA TN 4191 


55 



Figure 12.- Comparison of zero-lift wing-fus elage pitching-moment 

coefficients by two methods. 
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Figure 13 .- Comparison of flight and wind-tunnel zero-lift wing-fuselage pitching-moment 
coefficients. Vertical lines in circle symbols are values of gT; — from equation ( 17 ) . 
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(a) Center of pressure of wing aerodynamic 
load due to pitching acceleration. 





(b) Aerodynamic load due to pitching 
acceleration for both -wings. 


Figure 14.- Pitching-acceleration wing aerodynamic -load parameters as 
functions of flexibility parameter qm^. 
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Figure 15 . - Measured and calculated radii of gyration in pitch squared as function of 

flexibility parameter qm R /W. 
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